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ABSTRACT 

Rapidly rotating Be stars are observed as shell stars when the decretion disc is viewed 
edge on. Transitions between the two implies that the discs may be warped and pro- 
cessing. Type II X-ray outbursts are thought to occur when the warped disc interacts 
with the fast stellar wind. We suggest that tides from a misaligned companion neutron 
star can cause the observed effects. We make numerical models of a Be star decretion 
disc in which the spin of the Be star is misaligned with the orbital axis of a neutron 
star companion. Tidal torques from the neutron star truncate the disc at a radius small 
enough that the neutron star orbit does not intersect the disc unless the eccentricity 
or misalignment is very large. A magnetic torque from the Be star that is largest at 
the equator, where the rotation is fastest, is approximated by an inner boundary con- 
dition. There are large oscillations in the mass and inclination of the disc as it moves 
towards a steady state. These large variations may explain the observed changes from 
Be star to Be shell star and vice-versa and also the Type II X-ray outbursts. We find 
the tidal timescale on which the disc warps, precesses and reaches a steady state to be 
around a year up to a few hundred years. If present, the oscillations in mass and disc 
inclination occur on a fraction of this timescale depending on the orbital parameters 
of the binary. The timescales associated with the tidal torque for observed Be star 
binaries suggest that these effects are important in all but the longest period binaries. 
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1 INTRODUCTION 

Be sta rs are rapidly rotating close to their break up ve- 
locity (|Porterl [l996) . They are early type main-sequence 
stars which have shown Ha in emission at least once. They 
are variable in brightness and spectra which show broad 
Hel absorption and emission at either visual or UV wave- 
lengths. It is found that the emission, and so presum- 
ably their discs, are only temporary and so Be stars be- 
come B stars and vice versa. Be star discs form and dis- 
appear on a timescale of a few months to a few years 
ilTelting et al.lll998l: iBiorkmanl l200(il, |2002|) . More recently, 
iRaioelimanana. Charles fc Ud alski (201(]) observed the op- 
tical emission from Be X-ray binaries in the SMC and found 
superorbital quasi-periodic variations on timescales of about 
200 — 3000 d that are thought to be related to the formation 
and depletion of the circumstellar discs. 

We choose to model the Be star disc as a viscous de- 
cretion disc. Several other models have previously been pro- 
posed to describe the physical structure of the disc such as 
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the wind-compressed disc mo del and the magnetically wind- 
compressed disc model (see iPorter fc Rivinius! [2003. for a 
review). In the wind-compressed disc model the rotation of 
the star produces a wind towards the equator. The mate- 
rial from both sides collides at the equator and the shock 
produces a dense region confined to t he equatorial plane by 
the o pposing winds from each side (|Biorkman fc Casinellil 
Il993h . However, non-radial effects because of the rota- 
tion of the star ruled out this picture from dynamical ar- 
guments. It cannot reproduce th e observed IR excess or 
kinematic structure (|Porterl [1997), An attempt to revive 
this model was made with th e addition of magnetic fields 
(|Cassinelli fc MacGregoill2000D . However, numerical simula- 
tions of this model still do n ot explain the observations when 
stellar rotation is included (|Owocki fc ud Doula| [2003), The 
most likely alternative to the wind-compressed disc model 
is the viscous disc model that we have used in this work. 

The disc is modelled as a viscous decretion discs 
with the mass expelled from the neigh b ourhood of the 
Be s ta r itself dLee. Saio fc Osakil Il99ll : ICassinelli et al.l 

2002). iHanuschikl ~] 19961 ) finds that the d iscs are ro- 

tationally supported and lOuirrenbach et al (|l997l ) de- 
duced that the circumstellar envelopes around Be stars 
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are thin discs from interferomet ry of Ha emission. 
IWood. Biorkman fe Biorkmanl (|l997l ) showed that the discs 
have an opening angle of about 2.5° with polarimetry and 
this is consistent with a Keplerian disc. This is equivalent 
to an aspect ratio H/R w 0.04. 

Decretion discs have a different structure to accre- 
tion discs bec ause the inner boundary condition is different 
|Pringkj|l99ll ). In an accretion disc the inner boundary has 
zero torque and accretion is allowed on to the central object. 
In a decretion disc, the inner boundary prevents flow there 
on to the central object. The disc grows when the star exerts 
a positive torque on it or decays when there is no torque. 
If mass is injected on to an accretion disc, then outside of 
the injection radius, the disc is a decretion disc. Contin- 
ual addition of angular momentum allows the viscosity to 
transport angular momentum outwards though the disc and 
the disc mass moves outwards. There are mathematical disc 
solutions for the surface density of decretion discs which ei- 
ther have infinite mass or divergent radial velocity l|Okazakil 
2007). There is a third steady state decretion disc solution 
where there is no flow either in or out. Flow may be pre- 
vented in the outer regions of a disc by a torque, such a 
tides from a companion star and in the inner regions by a 
magnetic torque from the central star. 

If the disc is viewed edge on the Be star is seen as a 
shell star. Two Be stars, 7 Cas and 59 Cyg have shown two 
successive shell events, which have been explained by a disc 
that is tilted wi th respect to the equatorial plane of the star 
l|Humme]||l998h . The variation in emission line widths and 
profile shapes are then due to the precession of the disc. 
The idea that a disc might change its inclination to the line 
of sight is furthe r invest i gated through observations of 28 
Tau (Pleione) by Iffiratal |2007f ). This star also changes be- 
tween B star, Be star and shell star on a 35 yr timescale and 
the intrinsic polarization angle changes in phase with these 
variation s . The cause of the precession is also not clear but 
iHummell (| 19981 ) suggested that it might be induced by tides 
from a binary companion as most Be stars have a neutron 
star companion. 

Be/X-ray binaries show thr e e types of X-ray activ- 
ity l|Sterla. White fc Rosnerlll986l ; iNegueruela et ai1ll998h . 
Type I outbursts are the periodic X-ray outbursts that co- 
incide with periastron passage with a luminosity Lx = 
10 36 - 10 37 ergs 1 . Type II outbursts last weeks to months 
and show no orbital modulation but have a larger luminosity 
of Lx > 10 37 ergs -1 . It has been suggested that the Type II 
outbursts are linked to the interaction of the fast stellar win d 
with a warped disc l|Negueruela et al 1 l200ll ; lQkazakll200l| y 
The third type of X-ray activity is persistent low luminos- 
ity X-ray emission with Lx < 10 34 ergs" 1 . In this work we 
consider the warping and precession of Be star discs with a 
neutron star companion on an orbit that is misaligned with 
the spin of the Be star and attempt to explain this activity. 

In Section [2] we consider general properties of accretion 
and decretion disc structure. In Section [3] we consider prop- 
erties of a decretion disc model warped by a tidal torque 
and in Section [4] we make numerical models. In Section [5] 
we apply the models to observed Be and B star binaries and 
consider how they are affected by a tidal torque. 



2 DISC STRUCTURE 

We model the structure and evolution of warped accretion 
discs accordin g to the simplified formalism described by 
IPringlel l| 19921 ). We consider the disc to be made up of annuli 
of width dR and mass 2irT,RdR at radius R from the cen- 
tral Be star of mass M\ with surface density T,(R, t) at time 
t. The material is in Keplerian orbits with angular velocity 
Q, = \J GM\ I R? . The angular momentum at radius R is 

L = (GM!R) 1/2 U = LI, (1) 

where I = (l x , l y , lz) is a unit vector describing the direction 
of the angular momentum of a disc an n ulus w ith / = 1. 

We use equation (2.8) of [Pringld l| 19921 ) and add the 
binary tidal torque, Tyd, and the angular momentum source 
term, L acc , to give 
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The source of angular momentum by mass accretion on to 
the disc from the Be star is 

-Radd^add S(R — -Radd) jBe, (3) 



27T_R a dd 



where 



jBe = (sin 9, 0, cos 6) 



(4) 



is the direction of the spin of the Be star. We assume the 
material added has angular momentum aligned with the Be 
star spin and it is injected at a radius 7? a dd so that f2 a dd = 
^K-Radd)- We work in the frame where the binary orbital axis 
is (0,0, 1). We discuss the tidal torque in Section [3] 

There are two viscosities. First v\ corresponds to the 
azimuthal shear, the viscosity normally associated with ac- 
cretion discs, and secondly V2 corresponds to the vertical 
shear in the disc which smoothes out the twist. The sec- 
ond viscosity acts when the disc is non-planar. With the 
a-parametrisation the viscosities are given by 

Vi = etiC s H (5) 

l|Shakura fc Sunvaevl[l973T ). where i = 1, 2, c s = HQ, is the 
sound speed and H is the scale height of the disc. The di- 
mensionless par ameters ai ^ 1 an d Q2 m ust be determined 
experimentally. iLee. Saio fc Osakil ( 199ll ) suggest that with 
a — 0.1 it will take too long for the disc mass to build u p 
so we choose ai = 0.3 (e.g. I Jones. Sigut fc Porter] I2008T ) . 
iQgilvid (|l99Sl > found the viscosities to be related by 



U2 



CK2 

ai 



1 4(1 + 7a?) 



2a\ 4 + a\ (6) 

and more recently lLodato fc Pried (|2010l ) have verified this 
numerically. With ot\ — 0.3 then we find the corresponding 
02 = 2.66. We note that varying a affects only the timescales 
of the effects discussed in this work. 



2.1 Steady State Discs 

In this section let us consider steady state discs. We set 
dL/dt = then take the scalar product of equation ([2]) 



with I to find 



= 11 
RdR 



(7) 



because Tua-l = (see Section f3 . 1 [) . The tidal torque is zero 
when it acts in the same direction as the angular momentum 
of the disc at that radius. Equation has the solution 



u x E = A + BR~z. 



(8) 



where A and B are constants to be determined by the bound- 
ary conditions. The first term is the accretion disc term for 
which v\Y> — const (far from the inner radius) and the sec- 
ond term is the steady state decretion disc term which has 
llPringld fl99lh and zero mass accretion rate 
through the disc. 

The mass transfer rate through the disc is 



M = 2ttRE(-vr) 
where the radial velocity in the disc is 
3 d 



«R = 



f 1 \ rj/ 1 2 

[vxY.R^ -v 2 R\^\ 



dR 



(9) 



(10) 
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l|Pringlelll98"lh . The accretion rate is defined so that M > 
implies accretion, the material moves inwards, and M < 
implies decretion, the material moves outwards. The general 
steady state surface density in a disc is 



v M 
3tt 



1- 



(11) 



where E* and z/* are the surface density and viscosity at 
R — R± and the steady rate of mass transfer through the disc 
is M. The first term is the normal accretion disc term with 
z/E = at R — R* and so zero torque at the inner boundary. 
The second term satisfies «r = when 81/ dR is zero there. 
With the additional second term, we can vary the mass of the 
disc and the accretion rate through the disc independently 
even with a steady state disc. For a Be star disc that is 
truncated by the tides of a neutron star companion we take 
only the second term, the steady state decretion disc term 
(see Section [2.41 for more on the truncation of the disc). 



2.2 Surface Density 

The density of the disc is 



p(R,z) 



pa 



-R ] 6 ( H 



(12) 



where zq is the height of the disc about the equatorial plane 
and is a function of radius, R. For a flat disc Zo(R) = 
everywhere. By integrating the density over the disc height, 
z, we find the surface density of the disc to be 



pdz 



(13) 



with the constant 



E* = PoV^H*. (14) 

IWood. Biorkman fc Biorkmanl (1 19971 ) find that at the stellar 
radius the disc aspect ratio is H+/R+ ss 0.04. We choose a 
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10 -11 gcm -3 and n = 7/2 and we 



typical density of po 
find 



= 0.49 



H„/R t 



R* 



A, 



(15) 



g/cm 2 "'"V °- 04 J V101W Vl " n g cm " 2 
We justify the choice of n through the viscosity. 

2.3 Viscosity 

For an isothermal disc the sound speed c s = const and so 
the disc scale height is H = c s /fl oc R? . We use the a- 
prescription for the viscosities as in equation (0 and find 



Ri 



iRz. 



(16) 



where i = 1,2. With equation |(HJ we have viE oc 2 in 
a steady state decretion disc and with (|13[l we found E oc 
J? _n+ 2 . The viscosity must therefore be parametrised as 



(17) 



However, because Vi oc R? with the a-prescription, for a 
consistent model we must choose n = 7/2. We find the vis- 
cosities at the stellar radius to be 



= 1.69 x 10 



Mi 



2.66 



H+/R* 
0.04 



17M e 



R, 



10 R G 



(18) 



for i = l,2. We note that if the disc was not isothermal, then 
we would change the radius power law for the scale height, 
H, of the disc. The value of n decreases so that the surface 
density would fall off more slowly with radius. This does not 
affect the shape of the warped disc in the numerical models 
in Section [4] onl y the scale height and surface density of 
the disc. IPorterl (2003) found that with a smaller value of 
n he was able to better reproduce the observed continuum 
emission. 

The first viscous timescale describes the timescale over 
which mass flows through the disc and the second viscous 
timescale describes the how the disc smoothes out the warp 
and twist. The viscous timescales in the disc are 



Vi 



(19) 



and so we find 

33.2 r^r 1 

V2.66/ 

' Hj/Rj 
0.04 



10 R 



Mi 

17 Mr. 



(20) 



for i = 1,2. The viscous timescales increase with radius in 
the disc. For example, the minimum timescales at the inner 
edge of the disc at the stellar radius are t vi = 294 d and 
t„ 2 = 33.2 d. An accretion disc reaches steady state on the 
first viscous timescale. The disc smoothes out on the second 
viscous timescale. 
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2.4 Tidal Truncation of the Disc 



3.2 Orbital Parameters 



Providing that the orbital period is shorter than the vis- 
cous timescale, t vi , in the outer parts of the disc then the 
tidal torques from the companion neutron star can trun- 
cate the disc. In the circular orbit case, the tidal trunca- 
tion radius, R t , can be estimat ed as where ballistic orbits 
be gin to cross (Paczvriskilll977l). For example, in Figure 6 
of iMartin. Tout fc Pringld (|2008T ) this radius is shown as a 
function of mass ratio of the binary. For Be star systems, the 
mass ratio is in the approximate range M2/M1 = 0.1 — 0.5 so 
the truncation radius is R t < 0.5a, where a is the semi-major 
axis of the orbit. In the eccentric binary case, the disc ra- 
dius can be estimated by considering the effects of tidal reso- 
nances. As the eccentricity of the bin ary increases, the trun- 
catio n radius of the disc decreases (|Artvmowicz fc Lubowl 
1 19941 ). The neutron star orbit will be larger than the outer 
radius of the disc unless the eccentricity is close to 1 or there 
is a large misalignment angle. 

We note that this has implications for theories of 
Type II X-ray outbursts. It has previously been suggested 
that they occur when the Be star disc becomes sufficiently 
large that the neutron star passes through the decretion disc. 
This is only possible if the disc does not becomes truncated 
fast enough which may happen if the eccentricity is close 
to 1 or if the misalignment angle is large enough that the 
truncation is not efficient. lOkazakil |2007l ) suggested that 
the truncation becomes inefficient for large misalignments 
9 > 60°. In this work we consider smaller misalignments 
and note that observed misalignments of Be X-ray bin aries 
are relatively small (e.g. IMartin. Tout fc Pringlell2009h . 



Most Be star binaries have eccentric orbits and so the tidal 
torque will not be constant over an orbit. In an eccentric or- 
bit, the stars spend a short time close to periastron but here, 
the torque will be stronger than in a circular orbit. Because 
the orbital period of the system is much smaller than the 
timescale on which the tidal torque acts (see Section [5} we 
can average the torque over an orbit. In this Section we find, 
for the average torque over an eccentric orbit, the equivalent 
average binary separation. 

With Kepler's third law the binary semi-major axis is 



/ GMP 2 



(23) 



where P is the orbital period and M = Mi + M2 and so 
M \ 3 ( P 



a = 6.50 x 10" 



18. 4M 



24.3 d 



(24) 



From Kepler's first law, the separation of the binary stars 
in an eccentric orbit satisfies 



Rh 



q(l-e 2 ) 
1 + e cos <h 



(25) 



where e is the eccentricity of the orbit and <f> is the azimuthal 
angle around the orbit. With Kepler's second law 



(26) 



where cf> = d<j>/dt. We integrate this equation twice over an 
orbit to find 



3 TIDAL WARPING 

It is widely suggested that the warp and precession observed 
in Be star discs may be caused by a misalignment between 
the spin axis of the Be star and the orbit of the binary 
companion. Here we consider further the tidal torque. 



3.1 Tidal Torque 

The tidal torque in the frame of the binary orbit, of radius 
-Rb, is approximately 



Ttid = 



GM2RT, 
2RI 



J 3/2 



(e z .l)(e z x I) (21) 



l|Lubow fc Qgilviell2000l ; lQgilvie fc Dubusll200ll ) to first or- 
der in the angle of tilt between the disc and the orbital plane. 
Here e z is a unit vector in the direction of the binary orbital 
axis, M2 is the mass of the binary companion star and Rb is 
the binary separation. This tidal torque is averaged over a 
binary orbit. The Laplace coefficient of celestial mechanics 
can be approximated for small z by 



J 3/2 



15 2 175 4 , 

1 H z H z + . 

8 64 



(22) 



We note that it does indeed satisfy J.T«d = (see Section[2]) 
so that the tidal torque is zero when the disc and orbit are 
aligned. Thus even for large tilt angles the tidal torque term 
is likely to have this form, although with modified coeffi- 
cients. 



(27) 



The tidal torque found bv lOgilvie fc Dubusl l|200ll ) in equa- 
tion (|21[) has a factor of 1/Rt which we must replace with 
its average in time 



PJo R\ 



dt 



a 3 (l 

The average separation is thus 



(28) 



i? b = a(l-e 2 )5 (29) 

(see e.g. iHolman. Touma fc Tremainel fl997). In Section [5] 
we find the tidal timescale for observed Be star binaries and 
there we must use equation (|21[) replacing Rh with _Rb so 
that we do not underestimate the strength of the tidal torque 
in an eccentric system. 



3.3 Tidal Warp Radius 

The tidal torque balances the viscous torque in the disc at 
the tidal warp radius. Outside of the tidal warp radius the 
disc is dominated by the tidal torque and inside of this radius 
by the viscous torques. The tidal warp radius is 



Rti 



2v 2 *(GM 1 )2R b 
3GM 2 i??" 2 



(30) 
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l|Martin, Pringle fc Toutll2007l , [20091 ). Note that in deriving 
this radius, only the first term in the expansion in equa- 
tion (|22l) has been used. This is the same radius for accretion 
and decretion discs which have only a tidal torque because it 
is independent of surface density. We expect that well out- 
side this radius the binary torque term, coupled with the 
viscosity, flattens the disc and it aligns with the binary or- 
bital plane. If the disc does not extend up to the tidal radius, 
the torque can still have an effect on the disc and cause the 
disc to move towards alignment with the binary orbit even 
if they don't completely align. 

With the standard parameters used in the previous Sec- 
tion we find with n = 7/2 that 



R tid = 3.68 x 10 1 
Mi 



24.3 d 



/ Q2 \ ; 
V2.66/ 



17 M 



WR 



M 2 
1.4M S 

H^/R, 
0.04 



M 



18. 4M 



(31) 



If n — 4 (as we have used for our numerical models in Sec- 
tion |4]) then we find i? t id = 6.1 x 10 12 cm for these general 
parameters. With the small orbital period shown here the 
warp radius is within the disc. However for Be star systems 
of longer period, it may be outside of the disc. In this case 
the viscous torques in the disc always dominate the tidal 
torques. The disc is relatively flat, but may be tilted from 
the equator of the Be star and precessing. 

3.4 Tidal Timescale 

The timescale on which the tidal torque acts is equal to the 
viscous timescale at R — Rud by definition. It scales as 



itid = 



We parametrise this as 



R* ( Rtic 



£tid — ^tid* 



(32) 



(33) 



where the tidal timescale at R = R+ is 

. 3 

P 



itid* 



1.1(1-6 

Mi 



24.3 
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-)(—)' 
A V2.66/ 
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M 



18.4 M 



R* 



WRr. 



0.04 



yr. 



(34) 



The disc aligns with the binary orbital plane on this 
timescale. This is longer than the orbital period and so the 
disc inclination should eventually reach a steady state. 



4 EVOLUTION OF THE MISALIGNED DISC 

We solve equation J2l) numerically with a first order explicit 
method (|Pringld ll992) including the non-linear terms. We 
concentrate on the effects of the tidal torque, without which 
the disc would remain flat and aligned with the spin axis 



of the Be star. The size of Be star discs has been observed 
in a few cases. iGrundstrom fe Giesl l|2006h used results from 
a Ha monitoring campaign to find the Be stars that they 
considered had discs with an outer edge in the range 3 — 
Hi?*. In this work we choose the Be star to have an outer 
truncation radius of R ou t = WR*- We use 200 grid points 
distributed logarithmically between the inner radius Ri n = 
R* and the outer radius i? ou t- 

The inclination of the disc angular momentum to the 
binary orbital axis is 



i(R) = cos 1 (l z ) 



and the azimuthal angle is 



j{R) = tan" 1 ( f- 



(35) 



(36) 



The tidal precession can be parametrised as 

T m = n tid X L (37) 

fsee lMartin. Pringle fc Toutl l2009). where we choose coordi- 
nates so that 



fitid = 



t V2 [Rtid) \-Rtid 



R V 



(38) 



This term has the effect of precessing the disc about the e z 
axis when the disc in not aligned with the binary orbital 
plane. This corresponds to the first term in the expansion 
given in equation (|22|) . 

Initially the disc has very little mass that is flat and 
aligned with the spin of the Be star as given in equation Q, 
where 8 is the misalignment between the spin of the Be star 
and the binary orbital axis. Mass is injected on to the disc at 
a constant rate M at a radius of i? a dd = 1.1 R*- The mass is 
always added with an angular momentum aligned with the 
spin axis of the Be star which in all models is misaligned to 
the binary orbital axis by 9 = 25°. 

We need three boundary conditions at the inner and 
outer edges of the disc. We choose two conditions on the di- 
rection of the angular momentum at both edges, dl /dR = 
with \l\ = 1. We need one more condition on the surface 
density. In this work we consider three different surface den- 
sity inner boundary conditions that we describe in the Sec- 
tion [471] We use a zero radial velocity surface density outer 
boundary condition, «r, = 0, at R ou t because we assume that 
the outer radius is truncated by the tidal torque from the 
neutron star. Tid al torques truncate the edge of a disc over a 
small radius (e . g. Ilchikawa fc Osakilll994l : iMartin fe Lubowl 
l201lT l. lOkazakl (|2002l ) use a 3D SPH code to model a Be star 
disc that is coplanar with the neutron star orbit. They find 
that the disc is truncated if ai <C 1. With the neutron star 
companion, the dis c becomes denser than the disc around 
an isolated Be star. lOkazaki fe Havasakil (|2007l ) extend this 
model to a disc that is inclined to the binary orbital plane. 
They have material accreted on to the Be star disc for 50 or- 
bital periods and then they turn off the accretion. They find 
that there is little precession until the accretion is turned 
off. With a one-dimensional code we can evolve the disc for 
a longer period of time. 

We choose the first viscosity 



f! = 0.1 



f ) nit 



(39) 
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and the second viscosity 



where we have defined the time unit 



tv 2 (Ft* 



(40) 



(41) 



which with equation (|20p for our standard parameters t„, — 
58.8 d. With these viscosities the viscous timescales are con- 
stant in radius through the disc and the evolution proceeds 
at the same rate at all radii. This is slightly different from 
the prescription we found in Section [2] 



4.1 Inner Boundary Condition 

The inner boundary condition is an approximation to a 
torque applied close to the star that depends on the mag- 
netic field of the Be star. We assume that the magnetic axis 
is aligned with the spin of the Be star. We first consider a 
torque that is strong enough to prevent all accretion back 
on to the Be star. This is equivalent to vb, = at R t or else 



-( 

dR V / \R=R t 

l|Prmglelll99lh . Next we consider a disc with no additional 
torque from the star s o that the inn er boundary condition 
is E = at R = R* l|Pringlel Il98lh . Finally we choose a 
magnetic torque that varies with the inclination of the inner 
edge of the disc. The torque is strongest when the disc is 
aligned with the equator of the star but weakens as the inner 
disc moves away from the equator. We parametrise the inner 
boundary condition with 



: 



(42) 



d_ 

dR 



(uiER^ 



K 



cos(i, — 0) 



1 



(43) 



where = fl(Rf,), i* = i(R*), f2 s t ar is the angular veloc- 
ity of the Be star and K is a constant we can vary. We 
note that the zero radial velocity and zero surface density 
inner boundary conditions represent the two limiting cases 
of possible inner boundary conditions. We consider all these 
boundary conditions further in the Appendix and show that 
the evolution is the same with an extra torque term acting 
on the disc. 



4.2 Zero Radial Velocity Inner Boundary 

The zero radial velocity inner boundary condition is given 
in in equation (|42|l . In this case flow is prevented from both 
the inner and outer edge of the disc and so there is no mass 
loss from the disc. Mass is added to the disc at a constant 
rate and so the mass of the disc increases linearly in time 



Mt c 



(44) 



For example if the mass accretion rate on to the disc is 
M — 10 -9 Mq yr _1 and the second viscous timescale at the 
inner edge of the disc is t V2 (Ri,) = 58.8 d = t c . At a time of 
t = 243 t c the mass of the disc would be Md = 3.9x 10~ 8 Mq. 

In Fig [T] we show the inclination and azimuthal angle 
of the angular momentum vector and the surface density 
of such a disc which has a binary companion misaligned 



by 6 = 25° and a tidal warp radius of i? t id = 5 R* . Mass 
is injected at J? a dd and initially the mass moves outwards. 
Once it reaches the outer boundary where the disc becomes 
truncated it reaches a steady state where there is no flow in 
the disc. The surface density continues to increase steadily 
in time once it has reached its steady state of vYj oc R~? 
l|Prinridll9"9ir i. The tidal torque is strongest at the outside 
of the disc and so initially the outer parts of the disc align 
with the binary orbital axis and the inner parts follow. The 
tidal timescale at the inner radius is t t id(R*) = 11.1 tc (with 
equation I32p . At this time, the disc has become significantly 
misaligned with the spin of the Be star (left hand plot) and 
significantly twisted (middle plot). The surface density has 
reached its quasi steady state. The tidal timescale at the 
outer radius in the disc is t t id(10 R*) — 0.35 t c . At this time, 
the disc is still close to alignment with the equator of the 
Be star. 

In Fig. [2] we consider a disc with a larger tidal warp 
radius of i? t id = 20 7?*. This is outside of the outer edge of 
the disc. The viscous torques in the disc are always larger 
than the tidal torques. The surface density evolution is the 
same as in Fig. [1] because it is independent of the tidal warp 
radius. The disc is less warped and twisted with a larger 
tidal warp radius. It moves more like a solid body. The tidal 
timescale at the inner edge of the disc is tud(R*) = 89.4 t c . 
At this time the disc has moved significantly from the ini- 
tial alignment with the Be star equator. In the bottom right 
hand plot we show the evolution of the angular momentum 
vectors at the inner and outer edge of the disc in time. Ini- 
tially the disc has l y = and both vectors are the same, on 
the right hand side of the plot. The disc precesses about the 
binary orbital axis at the origin. 

From the two figures in this Section we conclude that 
the Be star disc aligns and precesses with the binary orbital 
plane on the tidal timescale at the inner edge of the disc. 
Once the disc mass is large enough, the additional torque on 
the disc from the injected matter becomes small compared 
with the torque from the binary companion in the larger 
mass disc so the disc always becomes aligned with the plane 
of the binary orbital axis, no matter what the warp radius 



4.3 Zero Torque Inner Boundary Condition 

The second inner boundary condition that we consider is a 
zero torque inner boundary condition normally associated 
with accretion discs 



= E* = 0. 



(45) 



As before material is continuously injected at i? a dd and it 
moves both outwards and inwards. In R > i? a dd the disc 
behaves as a decretion disc and in R < R ac \d as a normal 
accretion disc. 

In Fig. Owe plot the evolution of a system with _R t id = 
5 R+. In this example the disc mass reaches a steady state of 
Md = 0.36 Mt c . Because the disc mass is lower than in the 
previous case, the mass addition torque is relatively stronger 
compared with the tidal torque and so the inner parts of the 
disc remain close to alignment with the Be star equator while 
the outer parts tend towards alignment with the binary or- 
bital plane. We do not show the evolution of the angular 
momentum vector in this Section because there is very little 
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Figure 1. The disc has a warp radius of i? t id = 5 -R* and a misalignment of 9 = 25°. The inner boundary condition is tiR = 0. The 
evolution of the inclination of the disc angular momentum to the binary orbital axis (left), the azimuthal angle (middle) and the surface 
density (right) of the disc. The times shown are t = 0.037 (long dot dashed lines), 0.33 (short dot dashed lines), 3.0 (long dashed lines), 
9.0 (short dashed lines), 27.0 (dotted lines) and 243 t c (solid lines), where the time unit is defined in equation 1411 . The mass of the disc 
increases linearly in time. 



precession with the zero torque angular momentum inner 
boundary condition. 

4.4 Varying Radial Velocity at the Inner edge 

In this Section we consider the inner boundary condition 
that varies with the inclination of the inner disc to the equa- 
tor of the Be star given in equation (|43|l . Unless otherwise 
stated we consider the Be star to be rotating at break up 
velocity so that f2 a tar = fi*. In Fig. [4] we plot the inclina- 
tion, azimuthal angle and surface density of a disc with a 
tidal warp radius of Rtid = 5 R* and K — 100. Once the 
inner edge of the disc moves away from the Be star equa- 
tor, accretion on to the Be star begins. The disc reaches a 
steady state that similar to the one with the zero torque in- 
ner boundary condition (shown in Fig. [3]) but now is closer 
to alignment with the binary plane because there is more 
mass in the disc. 

In Fig. [5] we plot the evolution of the mass of the disc 
(left) and the inclination of the inner disc edge relative to the 
binary orbital plane (right) for discs with i? t id = 5, 20 and 
80 i?* and K = 10. There is competition between the tidal 
torque, the mass addition torque and the magnetic torque 
which causes oscillations in the mass and inclination of the 
inner disc. Eventually the disc always reaches a warped and 
twisted steady state. The larger the tidal warp radius, the 
larger in both amplitude and number the oscillations are in 
mass and inner inclination of the disc. With a larger warp 
radius, the steady state mass of the disc is larger and the 
inclination of the inner disc is closer to the Be star equator. 
The discs shown in Fig.[5]have warp radii of iZtid = 5, 20 and 
80 Ri, so their corresponding tidal timescales at the inner 
edge are 11.1, 89.4 and 715 t c . This is the timescale that the 
discs reach their steady state. The oscillations occur on a 
timescale shorter than this. 

Once the disc mass has built up to a steady state, the 
only way for the oscillations to reoccur is for the disc mass 
to be lost or the accretion on to the disc to slow or cease. For 
example, if the driving mechanism behind the mass ejection 
from the star is pulsation then it is possible that the oscil- 



lations will continue as the accretion on to the disc varies 
over time. 

In Fig. [6] we consider the effect of varying the constant 
K. We plot the surface density and inclination of the inner 
parts of the disc, and the total disc mass with varying K 
for Rtid = 20 R*. The larger K, the smaller the disc mass 
and surface density at the inner edge. The larger K, the 
larger the steady state inclination of the disc to the binary 
orbit. We also plot the evolution of the zero radial velocity 
inner boundary condition and see that when the tidal warp 
radius is outside of the disc then there are small oscillations 
in inclination of the disc as the disc approaches the binary 
orbital plane. Precession of the angular momentum vector 
is shown in the bottom right plot of Fig. [2] also. 

In Fig. |7| we plot the evolution of the angular momentum 
vector at the inner and outer edge of the disc for a disc with 
i?tid = 20 R+ and K = 10 which we can compare with the 
similar plot in Fig. [2] which instead had a zero radial velocity 
inner boundary condition. We see clearly the disc precesses 
about the steady state solution. 

In all of the models shown so far, the viscosity ra- 
dius power law index is n = 4 (see equation [T7|) . However 
IWaters. Cote fc Lamersl l|l987l ) and iDoughertv et al.1 1| 19941 ) 
matched the infrared excess of a large sample of Be stars and 
found the index to be in the range 2 ^ n ^ 5. An isotherma l 
disc theoretica lly requires n > 3.5 for outflow (|Porter|ll999l ) . 
More recently IJones. Sigut fc Porter! ([2007) used hydrody- 
namical simulations of the disc structure and predicted the 
range 3 ^ n ^ 3.5 in the inner regions of the disc. In Fig. [5] 
we show the effect of varying the viscosity parameters n 
and <xija\. Comparing the solid line (with n = 4), the long 
dashed line (with n = 7/2) and the dotted line (with n — 3), 
we find the smaller n is the longer the timescale of the os- 
cillations. However, qualitatively the results are similar. We 
note that when we vary n, we change the viscous timescales. 
In order to compare models with the same tidal warp radius 
but varying n, the second viscous timescale at the inner edge 
of the disc changes as t c (n) = t c (n — 4)(i?*/i? ti d) 4- ". Com- 
paring the long-dashed line (02/01 = 0.1) with the short- 
dashed line (q?2/ '<Xi = 0.02), we see that with a smaller 
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Figure 2. The zero radial velocity inner boundary condition for a disc with ittid = 20 R* and 9 = 25°. The evolution of the inclination 
of the disc angular momentum to the binary orbital axis (top left), the azimuthal angle (top right) and the surface density (bottom left) 
of the disc in time. The times shown are t = 3.0 (dot-dashed lines), 9.0 (long dashed lines), 27.0 (short dashed lines) 81 (dotted lines) 
and 243 i c (solid lines), where the time unit is defined in equation 1411) . The mass of the disc increases linearly in time. The bottom right 
plot shows the evolution of the angular momentum vector of the inner (solid line) and outer (dashed line) edges of the disc. At t = the 
disc is at l y = on the right hand side of the plot. 



value of olz/oli the oscillations become slightly larger and 
the steady state is closer to alignment with the binary orbit. 
However, the qualitative behaviour of the system is similar 
with varying 02/01. 



We have only considered Be stars that are rotating at 
break up with £7 s tar = fi*- We note that if fl sta r < fi* then 
the oscillation amplitude decreases and the disc moves to 
steady state more quickly. 



5 OBSERVED BE STAR BINARIES 

In Table [1] we show some Be and B star binaries that have 
been suggested to have a warped or tilted disc or have a mea- 
sured misalignment between the Be star spin and the binary 
orbital axis. There is one B star binary because it has a mea- 
sured misalignment and it is thought that B stars become 
Be stars and vice versa. We show the spectral type, the pe- 
riod, P, the eccentricity, e, and the measured inclination, 
6, between the spin of the Be star and the axis of rotation 
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Figure 3. The inner boundary condition is = for a disc with i? t id = and 8 = 25°. The evolution of the inclination (left), 

azimuthal angle (middle) and surface density (right) of the disc in time. The times shown are t = 0.037 (long dot dashed lines), 0.11 
(short dot dashed lines), 0.33 (long dashed lines), 1.0 (short dashed lines), 3.0 (dotted lines) and 81 t c (solid lines). The time unit is 
defined in equation ^411 . 
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Figure 4. The inner boundary condition is the one given in equation l|43|l that varies with inclination of the inner disc to to the Be 
star with K = 100. The disc has i? t id = 5i?* and 6 = 25°. The evolution of the inclination (left), azimuthal angle (middle) and surface 
density (right) of the disc at times t = 0.33 (dot dashed lines), 1.0 (long dashed lines), 3.0 (short dashed lines), 9.0 (dotted lines) and 
27.0 t c (solid lines) where the time unit is denned in equation 14H . 



Be/B Star Binary 


Spectral Type 


P/d 


e 


e/° 


4U0115+634 (V635 Cas) 


B0. 2 Ve 


24.3 


0.34 M 


unknown I 2 ' 


59 Cyg 


Ble I 3 ! 


28.2 


0.2 M 


unknown I 5 ' 


J0045-7319 


Bl V 


51.2 


0.808 I 6 ' 


25-41 m.M'l ] 


A0535+262 (V725 Tau) 


B0 Ill-Ve I 10 ' 


ill I 11 ! 


0.47 I 12 l 


unknown I 13 ' 


C Tau 


B2 III 


133 I 17 ' 




unknown 


0053+604 (7 Cas) 


B0.5 IVe 


203.59 


0.26 M 


w 25 I 8 ! 


28 Tau 


B8e 


218 


0.6 I 15 l 


59 I 16 l 


SS 2883 (B1259-63) 


B2e 


1236.7 


0.87 M 


> 55 PUS] 



Table 1. Binary Be stars and binary B stars that have either an observed misalignment or inferred warping. We show in col- 
umn 3 the orbital period, P, column 4 the eccentricity, e and column 5 the misalignment, 0, between the spin of the Be star 
and the bi nary orbital axis. Wher e t here is no measur e d mis ali gnment we give r eferences for the inferred warped , tilted or precess - 
in g disc. ^iRappaport et al.l Jl978h ' 2 ]NeguerueIa et al.1 j200 if) > 3 ]Harmanec et all 1120021) t 4 ]Rivinius fc Stefll i2000h l 5 lHummeJ l|l998|l 



t 6 jWex et all lll99Sl'l ^iHughes! jl999T) l 8 <Kaspi et al.l Jl996l> l 9 <Lai, BOdsten fc Kaspi 
|12 ]NeguerueIa et al l ||2000|) 113 fLarionov Lvutv fc Zaitseval j200ll ) l 14 lHarmanec et al 
l 17 lHarmaned (|l984l 'l l 18 ISchaefer et al.l l|2O10h 



19951) I 10 ISteele et al.l jigpsjl ^ilHutchinesj |l984 ) 
l|2000h l 15 lKatahira et al.l h99^ l l 16 lHiratal 120071 ) 
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Figure 5. Three models with the varying inner boundary condition in equation (1436 with K = 10. Left: The evolution of the total 
mass in the disc. Right: The evolution of the inclination of the inner edge of the disc. The spin of the Be star is at a misalignment of 
6 = 25° in all three cases. The dotted lines show a model with iJ t id = 5 R*i the dashed lines with i? t id = 20 -R* and the solid lines with 
Rtid = 80 R*. The time unit is defined in equation 141t . 
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Figure 6. Discs with i? t id = 20 R*. Left: The surface density at the inner edge of the disc. Middle: The inclination of the inner edge of 
the disc. Right: The total mass of the disc. In each plot the dotted line shows a disc with = (not shown in the left plot because it is 
zero). The solid line shows a disc with vn = inner boundary condition. The dashed lines show discs with the varying inner boundary 
condition with K = 10 (long-dashed lines) and K = 100 (short dashed lines). 



of the disc. Even if there is no observed misalignment for a 
Be star system, the formation mechanism of the companion 
neutron star by an asymmet ric supernova did most proba - 
bly produce a misalignment (|Martin. Tout fc Prin glc 2009). 
There is a large range of orbital periods and eccentricities 
in the Be/B star binaries. 

Ideally we should make an individual model of each 
system in order to draw firm conclusions about its evo- 
lution. However, in this Section we just consider the gen- 
eral properties of the Be star systems. We choose the 
mass of the Be star according to its spectral type (see 
lEggleton. Tout fc Fitchettl Il989h and then find the corre- 



sponding stellar radius from the mass (e.g. iTout et aljl996h . 
We list the properties in Table [2] 

In Table [3] we show the disc properties in each Be star 
binary. We find the semi-major axis of the orbit from equa- 
tion (124[1 and then the periastron separation is p = (1 — e)a. 
We find the average separation with equation (J29j) . The tidal 
warp radius is found with equation (|30|l with n = 4 and then 
the tidal timescale with equation (|32p . The tidal timescale 
tells us the timescale on which the disc reaches steady state, 
as shown in Section 14.41 Any oscillations of the disc mass 
and inner inclination of the disc occur on a timescale shorter 
than this. 
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Be/B Star Binary 


a/R* 


p/R* 


Rb/R* 


Rtid/ R* 


*tid»/yr 


4U01 15+634 (V635 Cas) 


9.3 


6.2 


8.8 


8.2 


2.1 


59 Cyg 


11.7 


9.3 


11.4 


12.7 


3.4 


J0045-7319 


17.4 


3.3 


10.2 


10.2 


3.3 


A0535+262 (V725 Tau) 


25.7 


13.6 


22.6 


54.4 


36.5 


0053+604 (7 Cas) 


38.5 


28.5 


37.2 


146.6 


161.4 


28 Tau 


67.1 


26.8 


53.6 


93.7 


44.9 


SS 2883 (B1259-63) 


165 


21.4 


81.2 


554.5 


825.2 



Table 3. The binary Be/B star binaries shown in Table [T] that have measured period and eccentricity. Here we compute the binary 
semi-major axis, a, the periastron separation, p, the average radius, R^,, the tidal warp radius, -Rtid and the tidal timescale at R*. We 
choose the mass and radius of the Be star system depending on its spectral type from Table [l] 




Figure 7. The evolution of the inner (solid line) and outer 
(dashed line) angular momenta of the disc with the varying in- 
ner radial velocity inner boundary condition with K = 10 and 
Rtid = 20 R+. At t = the disc begins with l v = on the right 
hand side of the plot. 



Figure 8. The inner disc inclination evolution for models with 
R tid = 40 .R* and K = 100. The solid line has a 2 /a ± =0.1 and 
n = 4. The long-dashed line has 02/^1 = 0.1 and n = 7/2. The 
short dashed line has 02/011 = 0.02 and n = 7/2. The dotted line 
has 0*2/01 = 0.1 and n = 3. 



The systems V635 Cas and 59 Cyg should have an evo- 
lution similar to that described by Fig. [4] as they both have 
small tidal warp radii. The mass and inner disc inclination 
evolution is similar to the dotted line in Fig. [5] They have 
tidal timescales of a few years and so should reach a steady 
state disc relatively quickly. The oscillations in mass and in- 
ner disc inclination should be minimal if the accretion rate 
on to the Be star decretion disc is constant. However, if the 
mass ejection rate from the Be star is not constant then the 
war ping and precession coul d continue over time. 

iNeeueruela et all (|200l[ ) observed the Be/X-ray tran- 
sient 4U0115+63/V635 Cas in optical, infrared and X-rays 
and that the disc warps, tilts and starts precessing. They 
found a quasi-period cycle with a period of 3 — 5 yr on which 
the Be star disc forms, grows and is lost. They observed two 
Type II outbursts before the disc became faint. The Ha line 
profile in this system is usually double peaked. However, 
just before the outbursts, the line became single peaked and 



this was sometimes followed by a shell event. They assume 
that the warping is because the disc becomes unstable to 
radiation driven warping. However, in such a short period 
binary with only a moderate eccentricity the tidal torques 
will certainly cause the disc to warp and precess if the sys- 
tem is misaligned. A s the warped disc prece sses, the line 
profile changes shape. iNegueruela et al.l (|200lT ) suggest that 
the warped parts of the disc will have a strong interaction 
with a fast stellar wind and stellar irradiation and so a large 
amount of gas falls on to the neutron star causing a Type II 
outburst and then the disc becomes faint because it has lost 
a lot of mass. We have found that the tidal timescale in 
the inner parts of the disc is about 2.1 years. This is the 
timescale on which the disc would reach steady state or else 
its maximum warping if the accretion on to the disc from 
the star is at a constant rate. However, with a varying mass 
ejection rate from the Be star, the disc will continue to warp 
and precess as the mass in the disc varies. 
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Be Star System Properties 


Be Star 










Spectral Type 


BO 


Bl 


B2 


B8 


Mass Mx/Mq 


17 


15 


12 


2.9 


Radius _R* / Rq 


10 


8.5 


7 


3.7 


Mass of Neutron Star, M2 
Disc Scale Height (H+/R+) 


1.4M 
0.04 









Table 2. Properties of the Be star systems we choose. We use the 
BO star parameters throughout the work and the other parameters 
for Table [31 



The B-star binary PSR J0045-7319 has a spin-orbit 
misalignment suggested by its orbital plan e precession 
l|Kaspi et al.lll99q ; lLai. Bildsten fc Kaspi 199^). In this case 
the B sta r rotates retrogradely with respect to the orbit 
i|Lailll996h . The system has a tidal timescale of 3.3 yr. This 
system should show oscillations in its mass and inner disc 
inclination on the timescale of a year or so if a disc was 
ejec ted. 

Haigh, Coe fc Fabregatl |2004h found that the disc in 
A0535+26 grows and decays on a quasi-cycle period of about 
1500 d by observing the change in quantised, IR excess flux 
states. The system has also shown type II outbursts in 1994 
and 2005 (|Coe et al.l[200r3 l. Because the system has only 
a moderate eccentricity, unless it has a very large misalign- 
ment angle, the disc will be truncated by tidal torques before 
it can reach the orbit of the neutron star. We find the tidal 
timescale in this system to be around 37 years and because 
the tidal warp radius is large, oscillations in mass and incli- 
nation should occur on a timescale of around the order of 
10 years. 

The Be star £ Tau has a companion that is yet to be 
detected. Because the eccentricity is unknown we do not 
compute the tidal timescales but note that the s ystem may 
be similar to A0535+26. Recentlv lSchaefer et al.l l 2010Ti used 
interferometry to observe the disc and suggested that the tilt 
of the disc is precessing. 

7 Cas showed two suc cessive shell phases over a period 
of 6 years (|Hummeli ri998). We find the tidal timescale at 
the inner edge of the disc to be 161 years. The tidal warp 
radius is very large, at 147 R t . This implies that the system 
will undergo large oscillations in disc mass and inner disc 
incli nation as the disc mass builds up. 

iHiratal l|2007f ) use polarimetry to observe the precession 
in 28 Tau (Pleione). Over the course of a hundred years, the 
star has changed from Be — > B, shell — > Be, shell — > Be with 
an activity cycle of 35 yr. It is currently in the Be phase. 
They find that over the time period 1974 to 2003 the disc 
axis precessed from 60° to 130°. With the Ha profile they 
suggest that the disc inclination has also changed drastically. 
They found the precession period of 28 Tau to be 80.5 yr 
but assumed a constant precession rate. Our model predicts 
that the precession perio d is 45 yr. The rigid body precession 
model of lLarwoodl (1 19981 ) found that the precession period is 
270 yr and 70 yr for disc radii of 10 J?* and 25 R* respectively, 
also in reasonable agreement. 

The tidal torque from the companion will have an ef- 
fect on all of the Be stars in listed in Table [3] apart from 
maybe SS2883 which has the longest orbital period and a 



very large eccentricity. The timescale on which the aver- 
aged tidal torque acts is very long. Because of the large 
eccentricity, a one armed spiral structure may fo rm at 
periastron passage (lHavasaki fc Okazakil [iooj . 120051 ). The 
fact that tidal torque s are negl i gible in this system has 
been n oted before by | Kochanekl (Il993l); iManchester et al.l 
l|l995l ); lLaT Bildsten fc Kaspil (|l995l ). IWex et al.l(|l998l ) sug- 
gests that the solution is a precessing orbit caused by the 
quadrupole moment of the tilted companion star. 



6 CONCLUSIONS 

We find that in Be star discs the torques from a neutron star 
companion are important and the discs can become warped 
and twisted if the system is misaligned (except in the sys- 
tems with the longest periods). Even though Be star's have 
eccentric orbits with their companions, the orbital period is 
short compared the timescale on which the tidal torque acts 
and so we can average the tidal torque over an orbital period. 
The discs are truncated at a radius smaller than the perias- 
tron separation of the binary providing the the eccentricity 
is not close to 1 and the misalignment is not too large. The 
size of the disc decreases with increasing eccentricity of the 
binary orbit. In most systems, it is unlikely that the neutron 
star will pass through the disc and so Type II outbursts are 
more likely linked to the warping of the disc. 

We have made numerical models of Be star discs that 
are assumed to be truncated by the companion and con- 
sidered solutions with three different inner boundary con- 
ditions that approximate the magnetic torque from the Be 
star. First we considered a magnetic torque strong enough 
to prevent all accretion back on to the Be star, or zero ra- 
dial velocity. In this case the disc mass increases linearly in 
time and the disc tends towards complete alignment with 
the binary orbital plane. Secondly we considered a disc with 
a zero torque inner boundary condition normally associated 
with accretion discs. The surface density at the inner edge 
of the disc is zero and the disc has a small mass. The disc 
moves towards a steady state solution with the inner parts 
close to alignment with the equator of the Be star and the 
outer parts become warped towards the binary plane but 
shows little precession. 

Finally we considered a magnetic torque that is 
strongest at the equator of the Be star (where it rotates 
fastest) and decreases as the disc moves away from there. 
The corresponding boundary condition varies the radial ve- 
locity at the inner edge of the disc with the inclination to 
the Be star equator. The mass addition and tidal torques os- 
cillate in dominance as the surface density varies. A warped 
and twisted steady state is reached once the oscillations die 
away on the tidal timescale at the inner edge of the disc 
which is of the order of a year up to a few hundred years 
(see Table |3I • On a timescale shorter than this the disc mass 
and inclination oscillate if the tidal warp radius is large. 
We suggest that these oscillations could be the cause of the 
shell events. The inclination and mass of the disc change 
on a timescale of a few, to a few tens of years until the 
disc reaches its steady state. If Type II X-ray outbursts are 
caused by a warped disc interacting with a fast stellar wind 
then the tidally warped disc solution could provide warping 
on a timescale similar to those observed. If the disc under- 
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goes periodic mass ejections followed by a decay of the disc 
mass, then the angular momentum distribution can be reset 
to being small and aligned with the Be star spin. As the disc 
mass builds up each time the oscillations will be observed. 

Observers should try to measure the misalignment be- 
tween the spin of the Be star and the binary orbital plane 
of more Be star systems. We note that the inclination of the 
disc is not necessarily aligned with either the binary orbit 
or the spin of the Be star and so the the two planes need to 
be measured independently of the disc inclination. 
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We can rearrange this to find 



/ = 



R- 



2(GMi 



vr + 







Ei?5 OR 



(A4) 



Al Zero Radial Velocity 

First we consider a torque that is strong enough to prevent 
all accretion back on to the star. This is equivalent to the 
inner boundary condition «r = at R = R+ no matter what 
the inclination of the inner parts of the disc relative to the 
Be star spin. In this case with equation ()A3|) we find the 
torque must be such that 

/= J— Af^EflO (A5) 

and because / falls off quickly with radius, this is approxi- 
mately equivalent to / = and 



= o. 



(A6) 



We can replace the source term in equation (|A1[) with this 
boundary condition. 

A2 Zero Torque 

Next we consider a disc with a zero torque inner boundary 
condition. We can choose / = and then the torque on the 
disc is the usual viscous torque 

T= -37n/iE(GMii2)i (A7) 

For zero torque at the inner edge we have E = at R — _R*. 



APPENDIX A: 
CONDITIONS 



INNER BOUNDARY 



Here we consider the torque on the inner parts of the disc 
from the rotating Be star. The tor que may be created b y 
the magnetic field of the star (e.g., iLivio fe Pringle! [l992]) . 
We need only consider a flat disc to investigate the inner 
surface density boundary condition because we have chosen 
the angular momentum to satisfy dl/dR = 0. The surface 
density evolution equation for a flat disc with a source of 
angular momentum is 



as 
dt 



RdR 
RdR 



3R? 



OR 
2AEi?i 



(Al) 



l|Pringlelll99ll ). where the input rate of angular momentum 
per unit mass is 



R J 



(A2) 



and f (R) ^ 1 measures the strength of the torque l|Pringlel 
Il99lf ). The function / is a strongly decreasin g function of 
radius dPapaloizou fc Linl Il984l ; IPringld Il99ll ). The radial 
velocity in the disc is 



vn 



Ri 







ERh dR 



(A3) 



A3 Varying Inner Radial Velocity 

The Be star rotates fastest at the equator. The third torque 
we consider varies with the inclination of the inner parts of 
the disc. The torque is strongest at the equator and falls 
off as the disc moves away from the equator. We choose a 
torque of the form 



f = C 





"-(<*- of] 


(f ) exp 


8 2 



(A8) 



where C and 8 are constant. We can vary the rate at which 
the torque falls off away from the equator of the Be star and 
the torque falls off quickly with radius from the star. 

If the constant, C, is large enough, this is equivalent to 
an inner boundary condition that takes «r = only at the 
equator, when the star rotates at break up, but allows for 
accretion on to the Be star as the inclination of the inner 
disc moves away from the equator of the Be star. At the 
equator = 8. As the inner disc inclination moves towards 
the pole, the viscous torque tends to zero and E* = and 
vr — > 00. Hence we can parametrise the radial velocity at 
the inner edge of the disc as 



t'R 



n*/n s - 



cos(u 



1 



(A9) 



where K is a dimensionless constant that we can vary. We 
cannot determine K without looking in more detail at the 
boundary layer between the disc and the star. The rotation 
rate of the Be star is Q s tnr and the Keplerian rotation rate 
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Figure Al. A disc with -Rtid = 20 -R* including the additional 
torque term given in equation IA8I and an inner boundary condi- 
tion £(i2*) = 0. We take the constant C = 15 and the upper solid 
line has S = 0.1 and the lower line has S = 0.4. The dotted line 
shows the case with the the inner boundary condition i>r = for 
comparison (the solid line in the middle plot of Fig. [6}. 



of the inner parts of the disc is S7* = f2(i?* 
condition at the inner edge of the disc is 



OR \R=R t 



R; 



n*/n s - 



cos(i* — 9) 



The boundary 



(A10) 



for cos(i* — 9) 7^ 0. If the star is rotating at break up then 

In Fig. I All we plot the inclination of the inner edge of 
the disc as a function of time for a disc with R t id = 20 J?* 
including the varying torque in equation (|A8fl and an inner 
boundary condition of E(i? t ) = 0. We see that this torque 
does indeed reproduce a similar behaviour to that with the 
varying inner boundary condition that we describe in Sec- 
tion The inner boundary condition is an approximation 
to the torque from a Be star that varies with height above 
the equator. However, the constants, C, 5 and K can only 
be determined by looking in more detail at the magnetic 
torque term. 



